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Abstract 



The Parisi-Sourlas mechanism is exhibited in pure Yang-Mills theory. Using the new scalar 
degrees of freedom derived from the non-linear gauge condition, we show that the non-perturbative 
0^ ' sector of Yang-Mills theory is equivalent to a 4D 0(1, 3) sigma model in a random field. We then show 

Cn ' that the leading term of this equivalent theory is invariant under supersymmetry transformations 

where x^ + 66 is unchanged. This leads to dimensional reduction proving the equivalence of the 
non-perturbative sector of Yang-Mills theory to a 2D 0(1, 3) sigma model. 

o, 

^ ■ 1 Introduction 

O 

.^-^ • There are now several new ideas on how to confine quarks inside hadrons in recent Hterature. The 

Q^ \ Seiberg-Witten paper IT] on the spontaneous breaking of N=2 supersymmetric Yang-Mills (SYM) to 

0^ ■ N=l SYM is supposed to implement the dual type II superconductivity model of confinement proposed 

by Mandelstam and 't Hooft 0. Last year, Fadeev and Niemi [bI derived an effective field theory for 

knothke solitons that is supposed to describe the infrared regime of pure Yang-Mills theory. In both 

f^ I approaches, confinement is essentially achieved through chromoelectric flux tubes, which provide the 

4j ' linear potential. These mechanisms break local Lorentz invariance because of the preferred direction 

defined by the flux tube. 

There is a third approach to confinement, which is based on the non- linear generalization of the 
is^ \ Coulomb gauge [Q . In a series of papers [g^ || 0| , this author argued the importance of generalizing the 

jj ■ Coulomb gauge to the gauge condition 

[d- D){d- A)^Q, (1) 

because field configurations which satisfy 9 ■ A" = /" 7^ 0, and [d ■ D)f = 0, are not gauge transformable 
to the Coulomb surface. In the non-linear regime (the linear sector being the Coulomb gauge) of the 
gauge condition ([|), the Yang-Mills potential can be written as 



Al - 7 ^^riS"' + £""7= + rhi-d^f" + ti). (2) 



(1 + 7^7) 
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As discussed in reference M , the effective dynamics of the new scalar field /" hints at non-pcrturbative 
effects. This was confirmed in reference H where the linear behavior of the instantaneous gluon prop- 
agator and the area law behavior of the Wilson loop were explicitly shown. The mechanism for these 
behaviors is not quantum mechanical but simply a statistical treatment of classical solutions of the ac- 
tion for pure /'s (all spherically symmetric /"(a;) in 4D Euclidean). This confinement mechanism then 
maintains local Lorentz invariance unlike the mechanisms that make use of chromoelectric flux tubes. 

Aside from not breaking Lorentz invariance, this mechanism relies on a straightforward, though non- 
trivial, difference between Abelian and non-Abelian theories. In an Abelian theory, gauge transformation 
is only a "translation" , thus all orbits pass through the Coulomb surface. Furthermore, the Coulomb 
gauge describes the physical degrees of freedom, the transverse photon. But in a non-Abelian theory, 
gauge transformation is a combination of "translation" and "rotation" that depend on the gauge field 
being transformed. This fact leads to complications like gauge copying phenomena [g| and the existence 
of gauge fields that cannot be gauge transformed to the Coulomb gauge |^, 0|. Hence, the proposal 
to generalize the Coulomb gauge to the non-linear gauge defined by equation |l|. Note that in the 
Abelian limit or in the short distance limit (weak coupling), the non-linear gauge reduces to the Coulomb 
gauge (/ — 0). But even before the coupling becomes very large, the fact that the non- linear term is 
proportional to p^ (p — momentum) while the linear term depends on p'^, the non- linear term would 
increasingly be more important in the long-distance regime. The non-zero / terms must now be taken 
into account and their dynamics provide the confinement mechanism. 

In this paper, we will prove a conjecture made in references H and W on the implementation of the 
Parisi-Sourlas mechanism |^ in pure Yang-Mills theory. The conjecture is based on two facts. First, 
the action for the pure / term is clearly non-perturbative because of (^) and it is infinitely non- linear. 

Second, the action for / is proportional to ((?/)■*, thus hinting that it can be written as (|j)^, where S 
is an action with a usual kinetic term (~ (9/)^). If indeed we can write the action for the / term as 
(Ij)'^, then its dynamics is stochastic, i.e., driven by a random field. Clearly this hints of the possible 
realization of the Parisi-Sourlas mechanism. 

The outline of this paper is as follows. In section ||, we give a detailed derivation of the path-integral 
in the non-linear gauge. Section ^ will prove the equivalence of the non-perturbative regime of Yang- 
Mills to an 0(1, 3) sigma model in random field. Section ^ presents the proof of dimensional reduction. 
The conclusion summarizes what the paper achieves. 

2 The path-integral in the non-Unear gauge 

The Yang-Mills path-integral in the gauge defined by equation |] is 

PI = hdA"^) 5{d ■ D{d ■ A)) detOe-^^''''^^\ (3) 

where O is the Fadeev-Popov operator 

O""^ = {D ■ dY\d ■ Df'^ - e'^^'^idf) ■ D"'^. (4) 

The important point about O is that it is a non-singular operator even though d ■ D is singular 0. 
The reason for this is that the gauge transform of a field configuration that satisfies d ■ A°- = f^ , and 
(9 • DY'^ f^ = 0, remains on the surface 9 • A" = /" but no longer in the horizon. This can be seen 
by considering an infinitesimal parameter A" — ef^, resulting in an A'^ = A'^ + eD"f'f^, thus satisfying 
d-A"" = f, but with {d-D''''{A')){d-A"') = -ee'^^'^D'lffdf.f ^ 0. This should rectify the inconsistent 
claim in reference (7) that the field configurations in the non-linear sector of the non-linear gauge never 
leave the horizon yet having a non-singular O. 

Let us introduce the scalar fields by inserting the following identity in equation g|: 

1= Jidnsir-d ■ A-). (5) 

The two delta functionals (see equations o and 0) imply that we can write 

llDff''^{A';,-d,^n-t% (6) 



where ft is a transverse vector field and ^o is a length scale introduced for dimensional reasons. Equation 
o can be solved for A't in terms of /" and i": 

Al = ^ . .^ {5'^' + e-'^ff^o'r + 9%rf){lld,f + t^). (7) 

Shifting t° +9^^/'' — » t° and rescaling g^o/" — ^ /"> '^^ find the expression given by equation ||. The 
new vector field t°, and the original potential A°, will now have divergence equal to -jif"'- 

What equation g means is that we traded the 12A° 's satisfying the three constraints given by equation 
y to the 12 t°j^s and 3 /"'s. The i°'s satisfy S^^JJ — -j^f"' (3 conditions). The extra degrees of freedom 
are removed by the constraints 

^a r abc i ahd rd fc acd rd rh i ra rd dbc 



1 + 7- / 



r (1 + 7 ■ 7)s'' - r(i + 7 • T)'^'^']^^/''*^ = o, 



Let us incorporate ft in the path-integral by inserting the identity: 



which are to guarantee d^Af, — -Trf^ 



M 



The path-integral now becomes 



PI = I {dtl){dn{dAl)6{d ■ D{d ■ A))6{d ■ A^ ~ -^n 



(8) 



«) s{t; - ti), (9) 

where ft is solved from equation (H). This means 



1 1 Sif^ - 1°) 

where 

.^=T'^Xj^ix-x'), (11) 

^"" = 7 ^— V ('5"'' + «°''r + /"/'')• (12) 

1 + 7- 7) 



1 / SA'^ \ 

X <5(A^ - ^'^"(-9^/'' + tl))det{0)det ^^j e'^-^^^l (13) 

Before integrating out A't, we can change the set of constraints from 

(j)" = {d- DY\d ■ A^) = 0, (14) 



x" = 5-^"-32r = 0' (15) 



1 
to 

r = p\ (16) 

1 



^-^d-e- -2 /^ (17) 



This is implemented by 



where 



5{r)Six'') = Si<jy"^)Six"')J-' 






J{^) = det 



gab j:ab 



5^{x-x')\ 



(18) 



(19) 



This Jacobian is 1, which can be shown by direct evaluation or by introducing fermionic coordinates 
and then doingthe integrations. Substituting the new constraints, evaluating the determinant given by 
equations ( |ll| , [l4 ), and integrating out A" we find 



PI 



1 



5'o 



(20) 



where it is understood that equation (||) is substituted in O and Sym- 

The Yang-Mills field strength in terms of /" and i° has the following form 



pa 



1 



M"^ „"M 



^1 + ^L + 5Q" (/;i), 



(21) 



where Z'ty, L't^, and gQ?i^{f;t) are zeroth-order linear and quadratic in t° , respectively. The action 
becomes 



^YM 



g 9 

2gLl^{f- t)Ql^{f- 1) + g^QlM; t)QUf^ *)] 



(22) 



3 Equivalence to a non-linear sigma model 

Consider the pure / path-integral 

PI{f) = [{dn det-\l + 7-7) detO e^^/'"'^^', 



(23) 



where O is the operator given by equation (m with f" = 0. To get the minus one power of det{l + f ■ f) , 

we changed the (1 + / ■ / ) factor of p" from minus two to plus one (see equation (||)). Why this change 
was made will be clarified later. In isolating the path-integral for pure / term (as given in equation (E3)) 
from the full path-integral, we are essentially claiming that the pure / dynamics is dominant and the 
remaining {f°",t°^) constrained dynamics can be treated as corrections. 
We will now show that this path-integral is equal to 



P/(/) = / {dD det 



srsf 



exp 



d^x 



dS 



where 



1 
2g 



s^^ I d^x v^'d^ra^f", 



^ab 



a £b 



r/ 



(24) 



(25) 



(26) 



(!+/•/) 
Equations (25 26) give the action for an 0(1, 3) sigma model in the non-liner form. Equation (e3) says 



that the non-pcrturbative regime of Yang-Mills theory is equivalent to an 0(1,3) sigma model in a 
random field. The rest of this section is devoted to the proof of this claim. 



The pure /° field strength can be written as 

z:,.^X""^dJ'd,r, (27) 



X"'" = _! , ( (1 + 27 • 7)6°'"= + 26"'' r - 26'"' f 

(1 + 7-7)2 ^ 

I n abd rd re __ ^^ acd rdb , bed ra fd\ 



(28) 



Equations (|27|j2^) give 
1 

4- '4^(1 + 7. 7) 

+ 6(1 + 27-7) [(a^r)'(/%/'')' - (df^rfd^f)^] } . 



z^ = {\) i ^ {(1 + 7 ■ 7r [{d,r d,rf (d.r o^rf 

* 1 + J • 7 r '^ 



(29) 



It is interesting to note that if we impose / -7 = 11 equation (|29| ) gives the interaction term of the 
Fadeev-Niemi action, which supposedly describes knotted strings 1^1 . This observation shows that indeed 
the /" dynamics describes the non-perturbative regime of Yang-Mills theory. 
The field equation of the non-linear sigma- model is: 

If we identify /" of above with the /" of Yang-Mills in the non-linear gauge, then we should have 

1 



°"'' ^ ^^f^lVf%f-nfV^^-^r%ff%f^] 



(31) 



linear term in i° 



Equation p^ comes from the non-linear gauge condition. Since we are considering pure / dynamics, we 
will put f^ — 0. Clearly we find that 

However, we can add a zero, a surface term, to the Yang-Mills action such that 

1 J d'xZ^ + J d^xd.H, -\jd'x {^^^ ' . (32) 

Inspection shows that 

a o £a\2 , Of i^ fa o fa\2 



+ l{d^rf{fd,f'f+6{d,,rfd^f'f (33) 

l2 



and equating corresponding terms in ( p2| ) yields 

2(7 • 7)' + 7(7 -7)2 -6(7- 7) -1 



a 



4(1 + 7-7)^ 
^ ^ 4(1 + 7-7)2' 



7 = 
5 = 

P = 



(7 -7)3 -11(7. 7) -7 

2(1 + 7-7)4 
(7-7)^ + 8(7-7) + 4 
2(1 + 7-7)4 
(7 • 7)' + 3(7- 7) + 1 



2(1 + 7-7)4 

To prove that the added term is indeed zero, we note that in 4D Euchdean space, 
giving 



d^x d,,Hn 



4., i^ V)^^ ,Q . ^1 i ^ .s. 



d'^y 



[x-yY 



(34) 



as the 4D Euchdean 



Since we will only be considering L^ fields, i.e., / d'^xA'tA't — finite, then A't ^ 

radius a; — > oo. From equation (0), we must have / ~ p^p^ as a; — > cxd. Shifting coordinates in equation 

(|34|), i.e., x^-yf, = y'^, we find 



(d^x{d-H) = i U d*y'yjL[{^d.H),-d^{d-H),y'^ 

J Jx^oo K^J y 



-do,dp{d-H)^y'^y'p- 



(35) 



rif^dS. 



We only have to look at the even number of integrands in y' because all the odd vanish. From the 
expression of {d ■ H) , we see that 



1 



1 



(d-H) 



Because of these behaviors, even though each y' integration diverges as y' — * cxo, the H factor that goes 
with it goes to zero faster as a; —> cxd resulting in a zero for each term. Thus, 



/ d'^x d^H^ = * 



H^ Hf^dS = 0. 



(36) 



To complete the proof of the equivalence to an 0(1,3) sigma model in random field, we need to 
establish the following determinant relation 



det'^il + J'-T) detO « det 



S^S 



(37) 



The proportionality constant may be infinite but it should be field / independent. At first glance, 
equation (p^) does not seem to make sense because O is a fourth-order differential operator while jt§t 
is only second-order. However, as we will show in the following arguments equation ( p7| ) is valid. 
It is important to note that we can write 



where 



S^S 



SfSf 
Mf = -2rd^f 



7j{a^ + Rif) + M^if)d^), 



1 + 7-7 



(38) 
(39) 



\2 + 7-7),rf'^-s'"' 



R = 5 -dafdar 

(i+r-r) 



(1 + 7-7 



1 + 



7-7)' 



(40) 



and rj is the "metric" given by equation (E6|). The factoring out of rj is important for two reasons. First, 
det 1] — det^^{l + f ■ f ) — det J-, and this cancels out the extra determinant at the LHS of equation 
(37). Second, it leaves n^ without a field-dependent coefficient, which in turn enables us to carry out 
the following identification: 

detO^deta^ det{n^ +M -d + R). (41) 

Equation (|4^) implies that the proportionality factor in equation ( p7| ) is det D^. It is hinted by the fact 
that for P = constant, O = (n^)^ and R = M^ = 0. 

Now we prove that equation ( [4l|) is true for any /°. Neglecting the t dependent term in equation 



(31), we find 9^tA° — 0. This docs not mean, however, that we are expanding about a field configuration 
on the Coulomb surface. Equation (||) clearly shows that this is not a case of background decomposition 
(where A = A + a) because t"^ is linked in a very nonlinear manner to /". Besides, we are not considering 
a particular field configuration /". Thus, if anything, the vanishing of the divergence of A° is purely a 
coincidence, but a welcome one because it leads to a significant amplification of O. 
Equations M) and (||), with i° = 0, give 



o = {ay + V^a^d^ + W^^d^d, + T^a^, 



(42) 



where 



yah 
rpab 



= -26"'"' Al,, 
= -5'^'A'^^Al 



A^K 



2e'^''^a^A^, 



= -^S''''Ald^A''^ + d^A^Al - e'^^'^u'^AI. 



(43) 
(44) 
(45) 



Equations (US 



5|) enable us to write 



where 



O^ [U^ + N -df, 



j\jab abc \c 



Since the determinant is invariant under a similarity transformation. 



(46) 
(47) 



det (D^ + AT ■ a) = det ( J"X° + N ■ d)J) ^det{U^ +Y -8 + P), 



(48) 



where 






Y^ = 2J-^d^J+J-^N^J, 



P = J-'U^J + J-'N^,d^,J. 



(49) 
(50) 



Using equations 



), equation mm is true if we can find a J such that 



tr ln{l +M -d + R) == trlnil 



n2 



{¥ -d + P) 



i_(r.a + p)ij(y.a + p)), 



(51) 





Figure 1: Equation (p2h in diagram form. 



where the trace is over space-time and isospin. Using the expansion for ln{l + x), we find that equation 
(15 1[) is equivalent to a set of infinite, global equations, 



tr 



1 

n2 



(Y-d + P) 



1 



1 

n2 



{M -8 + R) 



(52) 



with n = f , 2, . . . , oo, for the nine elements of J for each space-time point. Equation (B3) in diagram 
form can be written as in Fig || where the insertions on the left are M ■ d + P while those on the right 
are M ■ d + R. Because of the (1 -I- / • / ) factors in the denominators of the insertions, each term is 
essentially non-perturbative. 

Can J be determined by the infinite number of global equations given by equation (p3)? Naively, 
the answer is yes because we have infinite equations for the infinite unknowns (9 components at each 
space-time point). We can visualize the problem better if we discretize space-time (lattice formulation). 
The set of equations given by (p2[) becomes a set of non-linear algebraic equations for the values of 
J°'(x,). 

By choosing a suitable ansatz for J, we can always have the number of equations equal to the number 
of unknowns. It should be noted that both sides have divergencies arising from the same point limit of 
the Green function of D^. This must be handled carefully using a suitable regulator. In the continuum, 
what can be done is to solve for a J for ti = 1 and then verify the relationship for n — 2,3, . . .. For n=l, 
equation (|5^) can be solved by the local equations 



A" 
TDda 



2 ^ 

ryaa 



(53) 
(54) 



Using equations 



]), and the fact that iV°° = 0, we find 



J- 



-lab 



d„j'"' = -M° 



J 



-labr—,2 jba 



J 



-lab Ajbc 



1 



NTdnr"" = -R"". 



(55) 
(56) 



There are many possible solutions to equations ( p5| , p6| ). Assume the simplest ansatz 

r'^ai^S^' + Pif)/^/', (57) 

With / = ( / • /)2 the inverse exists under a very broad condition, i.e., c?[a + P f ■ j) ^ 0. The 
inverse is given by J-^^" = (i) 5°" - (|) ^j+^/'^Z". 

Equation ( pq ) leads to a first-order, non-linear, ordinary differential equation for a and /3 by equating 
coefficients of /"S^/", which is the only possible structure of M"" and Y^". Equation (pq), on the other 



hand, leads to two equations - a second-order, non-linear ordinary differential equation by equating 
coefficients of f^d^^f"" f^d^f^ and a first-order, non-linear ordinary differential equation by equating 
coefficients of dfj^f^dfj^f^. Seemingly, we have an overspecificd problem. Fortunately, the two first-order 
equations have the same structure and they only result in the following linear relationship between a 
and f3 

"(/) = /C(/)/3(/), (58) 

where 

^ ~fMf) ± VaHf)f' + SPif^ + 3/^ + 2)a(/) 

^ = Ml) ' ^''^ 

a{f) = 2/4 + 5/2-5. (60) 

We already see from equation lp^) the multiplicity of solutions even for the simpliest ansatz for J. 



Substituting (58) in the second-order equation, the problem simplifies tremendously, the differential 
equation becomes linear given by 

d^(3 



where 



-^+Bif)p + Cif)P = 0, (6f) 



4 f"^ IK' 2 f^AT' 1 
Bif) = -^ + ___ Uz t (Q2\ 

^■" 3/C + /2 /C SICiJC + p) /(1 + /2)' ^ ' 



C{f) 



f/C"(3/C + /2) JC'iSJC + p 



3 K[lC + p) 3/(l + /2)/C(/C + /2) 

(2/C -f- 2)f il-f)f 

3/C(f + /2)(/C + /2) 3(1 + /2) ■ ^ ^ 



Equations ( |6fy62|j63| ) are too complicated to have a closed form solution. However, a power series solution 
can be given because there are no poles in /. Remember that in the non- linear regime / 7^ and that 
/ is always real. This verifies the existence of a J that at least satisfies the one insertion condition of 
the trace In. Note that even at the one insertion level, the "equality" of the determinant is already 
non-perturbative. 

Finally, we answer the question, "Why not prove equation (H^) directly using A'^?" The answer is 
simple, X — thus we cannot have an insertion per insertion comparison of the tr Ins. The single 
insertion trace of F • 9 -|- P is contained in the higher insertion traces of N ■ d and comparison is difficult 
to make (but the terms are there but with different coefficients). 

4 Proof of dimensional reduction 

The path-integral given by (|23) can be written as 

PIif;w;T;r) = Jidn{dw'^){dT)idr)e-^, (64) 

where 

Because of the metric 7/°'', this action could not be derived from the supersymmetric version of S, i.e., 

A ^ fdWBde i77a6($)(a,.$'^a^$'' -i- 9e$''a^*'), (66) 

where 



Also because of 77, the supersymmetrized 5'($) is not invariant under the following transformation that 
leaves x'^ + 99 invariant, 



Xfj, — *• x^ 



9' = 9-2pe- X, 
9' = 9-2pe- X. 



(68) 
(69) 
(70) 



In equations ( |68| , |69| , [70| ) , e^,p,p are infinetisinial coordinate and grassman parameters. Dimensional 
reduction could not be proven using the non-linear action of 0(1,3). 

Actually, the paper of Parisi and Sourlas hints that we should use instead 



Sa^ J d^xC-d^ad^a - \d^,rd^r + A(fT2 -J -J - 1)}. 



(71) 



The equation of motion for the multiplier A(a;) implements the equivalence of Sa to S. We will show 
that the path-integral given by equation (p^) is equal to 



Piil.'yA) 



{dn{da){dX)Sl^-^)s^^^'' 



s\ 



Sa 



X det 



5'S. 



(72) 



where 



^6(l)iS(f)j 
(/°, a, A). The field derivatives are 



exp ■ 



^ / d^x 



6f- 



= a'r~2xr, 



^ = -nV + 2Aa, 
da 

f = "'-(1 + 7-7). 



It is straightforward to show that 






SS 



We will now give a "physicist's" proof of 

det ' 



Uiu^l 



det 



SScT SSo 



=0 



[srsfb 



(73) 
(74) 
(75) 

(76) 
(77) 



where the proportionality factor is anything as long as it is field-independent. Note that the matrix at 
the LHS of { \n\ ) is an "infinite 5 x 5" matrix while the RHS is an "infinite 3 x 3" matrix. However, 
direct inspection shows both determinants have leading terms (D^)^, thus equation ( |77| ) is not really 
surprising. Equation (76) can be written as 



5S 



SS„ 



Differentiating again with /°, we get 

6^S 



aXfb 



6f-6f 



S4>j 5f^ I SS„ ^SS„ ^Q 






The first matrix of the RHS of this equation is "3 x 5 (infinite)" , the second is "5 x 5" , while the last 
is "5 X 3" for consistency with the LHS which is "3 x 3" . Taking the determinant of both sides and 



10 



using the cyclic property of the determinant, we need to evaluate the determinant of the "5 x 5 infinite" 
matrix 



1 








■ 1 



fl 

(7 

sx 

■ W 



1 



fl 



sx 

w 





1 



fl 

(7 

sx 

w 




1 

fl 



sx 






1 

fl 

a 
SX 

W 



1 



1 



II 



1^ 



fl 

a 

fl 

a 

fl 

a 

7-7 



7-^ 






r_ 




a 


1 


fl 




a 


sx 




W 




sx 




7P 




sx 




W 




(7^ 

V 


5X 

s7 


<5A 


sx 


s7' 


s7 



sx 
W 
sx 
W 
sx 
W 



{S\x~x')) 



1 
Above we made use oi a — (l+/-T)') and A = ^ (i) D^ct — ^^9^/°9^/°. The important point 



is that by a series of manipulations, the fourth row or the fifth row can be made zero. Thus, 



det 



Sf-Sf 



= "0 • oo" det 



(78) 



The oo comes from the 5'^{x — x'). 

Effectively, we have established that the dynamics of f"" with action 5 in a random field is equivalent 
to a constrained stochastic dynamics of f^ , cr, A with action Sa- That this is true follows from the fact 



that we can always add ^ {j§^) and ^ (^") in the exponentials of (|7! 



'SS„ 



Next we complete the proof of dimensional reduction. First, we will exponentiate 5 (^-) by using 



{dip)S{ip ~ (pa)F{ip) 



C: 



{dip)exp{ 



8^P- 



1 



dx{(p - (pQf}F[ip) 



(79) 



+ ■ 



have exponentiated (S (^-) by 



It may seem that the use of the above approximation unnecessarily complicates things because we could 

ss„ 

So 

Um-^^,exA-\-Jd'x(^-^ 



Unfortunately, the presence of e invalidates SUSY except for e = 1. 

Second, we will exponentiate 6 (-rf-) by introducing an auxilary field Wx by 



J^)=j(dwx)eJ -- 



This was not done for ^- because a has a kinetic term (and A does not), which requires the presence 
of a ^w^ (will arise from dg^adg^cr)- 

Taking everything into account, we have 



P/(23) 



P/(| 
PUH 
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{dr){da){dX) det 



S'S, 



(80) 



X exp < — d X 



1 /6S„Y , 1 fSS^Y ^ 5Sx 



where the numbers in the PI refer to the equation numbers in this paper. Note that the last equation 
only involves an approximation, because we neglected the other terms in equation ( |79[ ) . Since the kinetic 
term of /" has the wrong sign, we momentarily shift if"" — » /", effectively rotating the 0(1, 3) symmetry 
to 0(4). Finally, we introduce the fermion fields for the determinant and the auxiliary fields Wa and Wa 
to complete the squares of the first two terms of the exponential to get 

P/(72) « j {dcp^){dw.,){d^,){d-^,) exp{-Ass}, (81) 

where 4>i = (/°,cr. A); Wi = {w'' ,w„,wx); ipi = {'iJj'',iJa,ip\); and 

/■ 4 f 1 2 1 2 < ^^<r , SSa 

Ass = dxl--w^--w,+w^—+w^ — 

(82) 

Equation ( p2[ ) can be derived from the supersymmetric version of Sa^ i.e., 
Ass^S,{^) = I d^xdMe\^d^<^„d^<^, + ]^de'^aC^'^a 



(83) 



where $° is given by equation ( |67| ) and 

$A = A + eijjx + i^x^ + ^^wx . 

Equation (|8^) explicitly shows invariance under equations ( |6q , |69| , |70| ) and dimensional reduction follows. 
Then we rotate back /° — > i/°, yielding again the 0(1, 3) sigma model in the linear form but in dimension 
reduced by two. 

This completes the proof that the non-perturbative regime of Yang-Mills theory in 4D is equivalent 
to a non-linear 0(1, 3) sigma model in 2D. 

5 Conclusion 

We have exhibited the Parisi-Sourlas mechanism in Yang-Mills theory. Since the starting point 
of the proof is the scalar field derived from the non-linear gauge condition, this paper proves further 
that the Coulomb gauge is not an appropriate gauge fixing in non-Abelian theories. Important field 
configurations will be missed in the path-integral as shown in reference B, where the linear potential 
was derived, and in this paper, where the Parisi-Sourlas mechanism was exhibited. 

It is also important to point out that the derivation of the linear potential is not quantum mechanical 
but merely a statistical treatment of a class of classical configurations, i.e., all spherically symmetric 
/°(x) in 4D Euclidean space. But as this paper showed, taking into account the full dynamics, including 
quantum effects, results in dimensional reduction and equivalence to an 0(1, 3) sigma model in 2D. 
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